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We give out a proposal of quantum simulation of pairing model on an NMR quantum computer. 
In our proposal, we choose an appropriate initial state which can be easily prepared in experiment. 
Making use of feature of NMR measure and the technology of the second (discrete) Fourier transfor- 
mation, our theoretical scheme can obtain the spectrum of paring model in principle. We concretely 
discuss the case in the concerned subspaces of pairing model and then, as an example, give out 
a simple initial state to get the gap of two the lowest energy levels in the given subspace. The 
quantum simulation to get more differences of energy levels is able to be discussed similarly. 

PACS numbers: 03.67.-a, 74.20.Fg, 76.60.-k 

Simulating a real physical system by a quantum computer (QC) was originally conjectured by Feymann Q. Later, 
this idea was confirmed by Llody in a two-state array and a general scheme for the quantum simulation was 
presented 0- 

NMR quantum computer is one of successful realizations of quantum computer so far Moreover, a four- level 

truncted oscillator |3( , a three-spin effective Hamiltonian ||j and the migration of excitation in an eight-state quantum 
system Q have been simulated. In recent, L.-A. Wu et al. reported an NMR experiment scheme perf orming a 
polynomial-time simulation of pairing models. Moreover, new works are coming forth continually 0, lid . In this 
letter, we give out a proposal of quantum simulation of pairing model on an NMR quantum computer. Two main 
features of our proposal are: (1) the choose of an appropriate initial state which can be easily prepared in experiment; 
(2) the using the second (discrete) Fourier transformation which can obtain the spectrum of paring model. We 
concretely discuss the case in the concerned subspaces of pairing model and then, as an example, give out a simple 
initial state to get the gap of two the lowest energy levels in the given subspace. The quantum simulation to get 
more differences of energy levels is able to be discussed similarly. In order to know what is the theoretical foundation 
of our scheme and whether our scheme can arrive at the needed precision, we also carry out some relevant research, 
sec our pre 

prints[nill2- 

In addition, after proposing our theoretical scheme we also finished its simple experimental 

implement |16|. 

Let us start from the spin-analogy of paring model Hamiltonian 0, 0] 



We first generally prove that if an appropriate working initial state is chosen rightly, the differences of energy levels 
of Hp can be obtained by using of two times Fourier transformations. 
Without loss of generality, a general working initial state can be written as 
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where = Y], bij\j) are the Hp's eigenvectors with the corresponding eigenvalues E l p , and \ j) are the standard spin 
basis. Here, we have used the fact that the eigenvectors |<^ mr ) of ff nmr with eigenvalues £^ mr are just \j) since H nmT 
in the laboratory is diagonal. 
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Considering the following evaluation 
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one obtains the density matrix 



P (t) = J212 a * bu ' a j b *n' ex P{-i ( E l - E l) T / h } 
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Recalling the procedure of the NMR measure, one can write the result, that is the NMR frequency spectrum, as 
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where cr^ = (a % x ± i<7*)/2, /?/ is the density matrix to be measured and ft means the Fourier transformation. It must 
be emphasized that the Fourier transformation is applied for the NMR measure time t. 
Now, substituting p(r) into the above equation, we have 



SnmrM « E E "'''""'A;' ex P{-i ( E p ~ H) T / h } S ( E Lr ~ ^ mr )Tr ( 1001 E 
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For an appearing peak in NMR frequency spectrum, for example E" mT —Eg ml , we can know its corresponding difference 
of energy levels (the values of a and (3) based on NMR Hamiltonian and the parameters of sample. Then, we carry 
out the second Fourier transformation and are able to obtain the -Hp's frequency spectrum 



Sp(E p ) cx ^2 aibi a a*b* j/3 6(E p - E J p ) 



(8) 
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In experiment, the second (discrete) Fourier transformation can be done from the dada that arc collected by measuring 
the areas (heights) of peaks E" mi — Eg mi at a series of evolution times Tj with H p . As soon as the E p — E 3 p are given 
out and the lowest energy level is also known, the spectrum of H p is just obtained. 

However, there are still three difficulties facing on us.(l) Which state is an appropriate working initial state for our 
purpose ; (2) Which E° ml - E^ ml . peaks can appear in NMR frequency spectrum for a given initial stste; (3) What i 
and j values correspond to the E p — E J p peaks in H p 's frequency spectrum. 

In order to solve them, we need some mathematical and physical preparations. 

Firstly, we should derive out the obvious expression Tr Ylk=l • Note the fact that the spin space 

can be divided into the different subspaces which correspond to the different numbers of spin-up states, that is 
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where the subspace n, i.e Sn" ' , is a subspace with n spin-up states |0). It 
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is clear that these subspaces contain the following basis 
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Further, based on the relations c + |0) = 0, er + |l) = |0), a |0) = |1) and a |1) = 0, we have that 
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In terms of these equations, it is easy to get 
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where in the last term n with possible values 2, 3, • • • , N — 1, i a ^ if, and i a = 1, 2, • • • , TV for any a = 1, 2, • • • , n. 
Secondly, we also need to analyze the structure of the eigenvectors of H p . From eq.Q) and the relations 

{o-i m) a { x l) +o-y n) o- { y ] ^ = +ic4 m) ) - > ( TO ^ 0, the y follow that |1) and |2 N ) must be Hp's two 
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eigenvectors respectively corresponding to the maximum and the minimum eigenvalues, which are denoted respec- 



tively by \v p ) and 



to Sn^ because that a + and er appear in pairs or do not appear in the various terms of H p . This implies that 



Moreover, if the arbitrary basis 



belongs to S^ N \ then H p 



also belongs 
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0, (m ^ n; m,n = 1, 2, • • ■ ,N — 1). Therefore 
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So we can denote the others eigenvectors of H p as v p ' \ — € with the corresponding 



eigenvalues -E^ and i^ n \j^ only take the sequence number of spin basis belonging to the subspace Sn\ for 
example i^(k) = 2 W - 2^"*), (fc = 1,2, • • • , AT). In other words, |i( n )), G S ( n N) . 

Now let us solve the difficulty one. If we only concern the spectrum E p ' and E p ' in the subspace 1 and 
subspace N — 1 respectively, we should choose such a working initial state that it does not include any 
2, 3, • • -N — 2). Obviously, the appearing NMR frequency spectrum will be simplified as at most 2N peaks 
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where in the summation ^ , i, j only take over 1, 2( Ar ~ fc ) + 1, 2 N - 2 ( - N ~ k \ 2 W ; (fc = 1, 2, • • • , TV). In this simplified 

case, the NMR spectrum only include the differences of energy levels in the subspaces 0, 1, N — l and N. In particular, 
if the working initial state is taken as c 



then E l p ' - E% can be obtained from the second Fourie 



transformation unless no any peaks appears in the NMR frequency spectrum. 

Then, let us solve the difficulty two, that is, how to guarantee the NMR frequency spectrum has those needed peaks. 
The key matter is that we have to take an appropriate working initial state including the kets |1), |2( JV— fc ) + l)(e 

S$_ x) ) or(/and) \2 N - 2^- k )){e S{ N) ), |2 W ), that is |^ ni ) = fli|l> a^-^vf^) or l^ini) = «2« \2 N ) + 

°i< 1 ) v p ^ Respectively, for the two cases we have 

N 



5*nmr(w)cx^ a^N-Db^N-D^N-^^alb^expi-^E^ E^j t}5(E^ t - i£ mr ) (14) 
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Since a i (N-i)6 i (w-i) j2 (w-fc)+ia*6*i or a 2 iv& 2 N ;2 iva* (1) &* (1) 2 «-2< N - fc ) are a ^ n0 * zer0 ! NMR spectrum must have some 
peaks appearing. 



The last difficulty is easy to be solved. In fact, setting IV'ini) = a2 N \% N ) + a 2 « _ 2 (jv-k) 
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E p from H p s spectrum and then obtain E p ~ in terms of E p = — i_, m=1 e m /2. If setting 



, we can read 
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(fci 7^ fc 2 ), we can read E p 



2 N_ 2 (N-k 1 ) 
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and Ef ~ 22 - <N - fe2) - , and then obtain tff -^"^ _ E f (get ^- 2 C-».) > £2«-2<— > if fc2 > fa) 
Therefore, in principle, we always can find out the energy levels and their differences in the subspace 1 of H p . Likewise, 
one ought to able to obtain the other energy levels in the other subspaces by the different chooses of the working 
initial states. 

It must be emphasized that it is interesting what is the physical meaning of E p — E p m theory 

and what is an appropriate working initial state to obtain E p — E p m experiment. 

Actually, by the submatrix diagonalization of spin-analogy of pairing model 1 1 1 1 ] and the numerical calculation |12| , 
we have found that the relation when N is large enough 



(£ 2 V 2 + A 2 ) 1/2 - feV: + A 2 ) 1 / 2 « £f - 2 - Ef-' 
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where E p 2 — E p 1 is the difference of two the lowest energy levels in the subspace 1 of H p , A is the solution of 
the energy gap equation 0, 0] 



while £ m comes from 
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Here 7^, 7 m are the quasiparticle creation and annihilation operators; e s is the ground state energy of superconduc ting 
system. Hamiltonian (|18fl was obtained by use of Bogoliubov transformation and mean field approximation ^[ Il5| 
to the BCS Hamiltonian iJecs = Sm=i( £m — £ F)(n m + n_ m )/2 — V^2m ;=i c tn c -m c -l c l- This result is consist with 
the known conclusion. Thus, we know the physical meaning of E p _2< x> 



2 N_ 2 (lV-fc 2 ) . 

-E; m theory. 



In order to see what working initial state should be chosen to obtain _E 2 
us consider an example when k 2 = 1, k\ = 2. It is easy to prove that 
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in experiment, let 

(19) 
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where we have defined 
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Here, i =/= j and i, j = 1, 2, • • • , AT. Obviously, are not completely independent each other. In practice, we can fix 
i = l and j = 2, 3, • • • ,N and obtain (AT — 1) linearly independent Note that when all e m are equal, so-called 

anti-W state \ W) corresponds to the lowest energy level in the subspace 1. When e m are different for the different 
m and — is a very small positive parameter, we know, from the perturbation theory, that the lowest energy 
level is still related with \ W), and the second lowest energy ought to relate with \ui2)- Thus, in order to make the 



working initial state to contain Vp 



(2 J 



and v, 



,( 2 ' 



(0) 
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definitely, we should take it as 



\2 N ) + \W) + \u 12 )) 



(22) 



Thus, the final -Hp's frequency spectrum includes consequently the peaks 



E p 



and R 2 "- 2 - 



Obviously, 



such an initial state is simple and easy to prepare in NMR experiment. 

In the summary, our scheme of quantum simulation of paring model on an NMR quantum computer has four 
steps. (1) prepare the initial state based on the quantum simulation purpose, for example 1)22(1 in order to obtain 
Ep N ~ 2 — H 2 ™ -1 ; (2) use NMR pulse series to implement evolution e~ lH P T / h 0,^J at a series of times t^; (3) carry 
out NMR measure (the first Fourier transformation) to clollect data which are a series of amplitudes corresponding 
to the evolution time Ti by H p ; (4) make the second (discrete) Fourier transformation for the collecting data and 
then obtain -Hp's frequency spectrum. Newly, we have finished the experiment of quantum simulation for the simplest 
system of two qubits. More experiment detail is put in our another paper 

Comparing with the known scheme [8j and our one, obviously, the working initial state in the former are prepared 
by the process of a quasi-adiabatically evolution, but our scheme does not need such a process and chose directly an 
appropriate working initial state which can be easily prepared in experiment. In terms of feature of NMR measure 
and the second (discrete) Fourier transformation, our scheme can obtain the spectrum of paring model in principle. 
It is different from the scheme in ref.@ where the second Fourier transformation was not used. In addition, because 
we chose an appropriate working state, the rotation step is not needed. However, it must point out that for every 
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time evolution, we use the same serial of pulses as the ref. |8( . These features of our proposal leads that the simulating 
paring model on an NMR quantum computer is actually feasible and really complete. Moreover, we concretely discuss 
the case in the subspaces 1 and N — 1 of pairing model and then, as an example, give out a simple initial state to get 
the gap of two the lowest energy levels in the subspace 1. The quantum simulation to get more differences of energy 
levels is able to be discussed similarly. This working is on progressing. 
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